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ABSTRACT 

Since HCN emission has been shown to be a hnear tracer of ongoing star formation 
activity, we have searched for HCN [J — 1 -^0) emission from known GMCs in 
the nearby galaxy M33. No significant HCN emission has been found along any of 
the lines of sight. We find two lines of sight where CO-to-HCN integrated intensity 
ratios up to 280, nearly a factor of 6 above what is found in comparable regions of 
the Milky Way. Star formation tracers suggest that the HCN-to-star formation rate 
ratio (Lhcn/-^*) is a factor of six lower than what is observed in the Milky Way (on 
average) and local extragalactic systems. Simple chemical models accounting for the 
sub-solar N/0 ratio suggest that depletion cannot account for the high CO-to-HCN 
ratios. Given HCN formation requires high extinction {Ay > 4), low metallicity may 
yield reduced dust shielding and thus a high CO/HCN ratio. The turbulence and 
structure of GMCs in M33 are comparable to those found in other systems, so the 
differences are unlikely to result from different GMC properties. Since lower CO-to- 
HCN ratios are associated with the highest rates of star formation, we attribute the 
deficits in part to evolutionary effects within GMCs. 



INTRODUCTION 



jSFR 



All local star formation is invariably associated with molec- 
ular gas, but recent studies have clarified that star formation 
is associated exclusively with de nse (n > 10^ cm~^) molec- 
ular gas from the s cales of cores (lLadalll992l:IWu et al.ll2005l . 
l20ld ) to galaxies (|Gao fc SolomonI 12004 iGao et al.ll2007^ . 
Most of the gas typically observed in studies of molecular 
clouds on large scales (> 10 pc) is found at lower densi- 
ties (n ~ 10^ cm~'^), but dense gas studies suggest that the 
low-density gas is not directly involved in the star formation 
process. This clarification is particularly important for stud- 
ies of the galaxy-scale star formation since the molecular gas 
content of a external galaxy is usually traced with emission 
from low-J transitions of ^^CO, which arc excited in the 
low density molecular gas which traces the bu lk of material 
found in molecular clouds (e.g.. IPineda et al.ll2008l '). 

Relatively few large-scale studies of the dense molecular 
gas content of galajdes have been conducted. The majority 
of mass in the dense gas is found as molecular hydrogen, 
which emits no radiation at the densities and temperatures 
typical of star forming clouds. Observations rely on emission 
from dense gas tracers, which are excited to emission only 
at the high densities. Typical tracers include line emission 
from HCN, N2H+, NH3 and HCO+, aU of which have critical 
densities of order 10^ cm~^. Whil e the relative merits of 
these tracers have been compared l|ShirIev et al ][2008). aU 
seem to produce similar results: the dense gas that they trace 
is linearly related to the star formation rate 



M, 



pc- 



(1) 



These results even seem to hold for high-redshift systems 
(Gao ct al. 200^. 

Even though low-density molecular gas does not par- 
ticipate directly in star formation, there has been exten- 
sive work demonstrating that low-density gas is indeed re- 
lated to star formation on a galactic s cale. While originally 
demonstrated for galaxies as a whole (|Kennicuttt. 1998) . re- 
cent work has shown the relationship holds for molecular 
gas o n scales of > 500 pc jW ong fc Blitz| l2002l : iHever et all 
2003: lKennicuttll2007l : bigieTe t al. 2001)7 These studies find 
SspR oc Sh2i where a = 1.0 to 1.4 depending on analysis 
methods used. The apparent difference in star formation re- 
lationships for low- and high-density molecular gas can likely 
be attri buted to cha nging dense-gas fractions in a galaxy's 
ISM {G ap fc Solomo n 2004). Some authors argue that the 
index of the star formation law can be predicted from the 
average density of the med ium compared to the critical 
density of tracer in que stion (|Krumholz fc Thompsorj|2007l : 
iNaravanan et al ][2Q0i), with additional considerations re- 
garding excitation of higher levels. Even with these consid- 
erations, the (1 0) line of HCN should provide a linear 
tracer of dense gas mass and an index of unity on the star 
formation law. Thus, a compelling model emerges where the 
star formation process in a galaxy can be decoupled into two 
steps: the formation of dense gas within a molecular cloud 
and the subsequent formation of stars in that high-density 
gas at a constant rate (or a constant fraction per free-fall 
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time, iKrumholz fc Thompso"nll2007l '). The apparent changes 
in the star formation efficiency (Eh2/Sspr) can then be 
attributed to change in the molecular cloud properties in 
galaxies, specifically the fraction of the cloud found at high 
density. 

Comparing the CO and IfCN emission from individ- 
ual molecular clouds can test whether there is significant 
change in the fraction of dense gas within individual molec- 
ular clouds. Some studies have been undertake n for Milky 
Way molecular clouds (e.g., 'Heifer & BlitzJ'l997l), suggesting 
that dense gas fraction of molecular gas does indeed change 
as a func tion of distance f rom t he centre of the Galaxy. The 
study of iBrouillet et all (|2005l ') examined the CO-to-HCN 
intensity ratio across the disk of the massive galaxy M31, 
finding a radial decline in the ratio with a maximum ratio 
of Ico/^HCN = 125. It remains unclear based on these re- 
sults the degree to which the changes in the CO-to-HCN 
ratio result from changes in molecular clouds properties or 
whether the relationship between the tracers and the under- 
lying H2 is changing. 

The nearby galaxy M33 has been not ed to exhibit a par- 
ticularly high star formation efficiency |Hever et al ] I2OO4I : 
iRosolowskv et al.|[2007l : ICardan et al.ll2007l ') based on obser- 
vations of CO emission. In the context of the above star 
formation model, this would imply that the dense gas frac- 
tions of the molecular clouds in M33 would be higher than 
that of the typical galaxy. Other circumstances may ex- 
plain the high star formation efficiency of the galaxy since 
M33 is distinct from M31 and the Milky Way in two ways. 
First, it has a significantly lower m etallicty outside the 
central region (12 -I- log(0/H) = 8.4 iMagrini et all [2OO7I : 
IRosolowskv fc SimonI I2OO8): and second, it has a signifi- 
cantly lower mass (Mmss ~ A/mw/IO). Given its proxim- 
ity and previous study it is possible to assess how tracers of 
high- and low-density molecular gas relate on the scale on in- 
dividual giant molecular clouds. To this end, we have under- 
taken observations of CO and HCN toward giant molecular 
clouds (GMCs) in M33 using the IRAM 30-m telescope. 



2 MOLECULAR LINE OBSERVATIONS 

We observed the positions of four Giant Molecular Clouds 
(GMCs) in M33 using the IRAM 30-m telescope located in 
Pico Veleta, Spain. The observations were carried out from 
16th to 20th of July 2007. HCN (1 0) and "CO (2 1) 
were observed simultaneously with the VESPA back-end in 
the 80 kHz resolution mode. The same configuration was 
used to observe ^^CO (1 -5> 0) and ^^CO (2 -s- 1) simultane- 
ously. All the observations were taken in position-switching 
mode. Weather conditions were good for summer observa- 
tions and pointing observations imply a pointing accuracy 
of < 5". Of note, this accuracy is small compared to the 
89 GHz beam (HCN(1 0)) but is more significant when 
compared to the 230 GHz beam ("CO(2 1)). Opacity 
and system temperature for each of the transitions is given 
in Table [T] and integration time {tint) and per-channel rms 
(cr) for the HCN/^^CO configuration is given in Table [1 

Conversion from the measured TX to Tmb is carried 
out as Tjnb = TX Fgf f / Beff ■ The forward {F^ff) and beam 
(B^ff) efficiencies were determined by measurements made 
by observatory staff. A summary of the observation's pa- 



Table 1. Observational parameters 



Line 




Preq. 






fpWHM 


/rp \ 

\-t sys/ 








(GHz) 






(arcsec) 


(K) 




1^00(1 -> 


0) 


115.2712 


0.95 


0.75 


22 


250 


0.28 


12C0(2 ^ 


1) 


230.538 


0.91 


0.52 


11 


250 


0.11 


"CO(2 ^ 


1) 


220.3986 


0.91 


0.57 


12 


250 


0.11 


HCN(1 -5- 


0) 


88.63160 


0.95 


0.78 


28 


110 


0.05 



rameters is presented in Table [T] and the spectra are shown 
in Figure[2] The data were reduced using the CLASS90 soft- 
ware. 

We selected four targets from the CMC catalog of 
IRosolowskv et all l|2007l 'l based on strong CO emission with 
a range of star formation properties. We summarize the ob- 
servations of the GMCs in Table [5] and their locations are 
shown in Figure [1] We did not detect HCN (1 — > 0) emis- 
sion in at least two of the four lines of sight. M33GMC 
1, the pointing closest to the center of the galaxy, shows 
a modest detection at a low level, as does M33GMC 76. 
For each line-of-sight, we measure the integrated intensity 
/hcn = X/ FmhSV where 5V is the channel width, where 
the sum is carried out over the expected velocity range of 
the HCN line assuming it has the same intrinsic line width 
as the CO line but accounting for hyperfine structure. We 
establish a 3a upper limit as 



/hcn < 3(7hcnAV [(1 - N/Ntot)] 



(2) 



where (Jhcn is the rms noise fluctuation in the spectrum, 
AV is the velocity width of the HCN line, A'" is the number 
of channels across AV and Ntot is the number of channels 
in the entire spectrum ([Matthews fc Gcioll200ll ). 

We further process these observations to place them 
on a common standard and supplement the results with a 
significant amount of ancillary data, the results of which 
are given in Tabled The integrated intensities are given as 
7io, hi and /i3 for the ^^CO{l -> 0), ^^00(2 ^ 1), and 
^^CO(l 0) lines respectively. For comparison to other 
work, molecular line intensities are scaled to luminosities for 
the lines by multiplying by the projected area of the beam 
for that tracer. 

For a rigourous comparison of line intensities, all data 
must be sampled with the same resolution. To properly es- 
timate the CO-to-HCN ratio, we convol ved the fully sam- 
pled, high resolution ^^CO(l —>■ 0) map of lRosolowskv et al.l 
( 2007, ) to the 28" resolution of the IRAM beam at 89 GHz 
and sampled the map at the observed positions. These val- 
ues are reported as 7io,ro7 and /io,convoi in Table [2] for pre- 
and post-convolution respectively. We checked the relative 
calibrations of the map by convolving the same map to the 
22" resolution of the IRAM ^^CO (1 0) observations and 
comparing the observed and predicted integrated intensities 
finding them consistent to the 10% level expected from the 
amplitude calibrations of the two maps. 

We supplement the molecular line data with infrared 
and optical band observations which are indicators of the 
star formation rate. To derive the star formation rate, we use 
combined Hq-|-24 /x m emission derived from data from the 
Local Group Survey (|Massev et al ]l2007l) and Spitzer/MIPS 
observations of the galaxy ( Polomski et al.l |2006| ) respec- 
tively. A continuum subtracted Ha image was generated 
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Figure 1. Locations of tlie four IRAM lines of sight overlaid on 
an integrated intensity of CO (1 — > 0) emission. Three contours 
of constant galactocentric radius arc indicated with ellipses rep- 
resenting -Rjipc = Ij 2 and 3 kpc. 
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by subtracting off an optimally scaled R band image from 
the Hq filter and s caling to photometric units based on 
iMassev et al] 1120071'). The M I PS da ta w ere from the cali- 
brated images of Hinz et al] l|2004 ) and iTabatabaei et al] 
l|2007l ). Each image was convolved to the 28" IRAM beam 
and sampled at the locations of the IRAM pointings. The 
lo cal star formation was established using the prescription 
of iKennicutt et al.l |20o3): 

{Mq yr"^) = 7.9 x 10"*^ [Luc + 0.039 U2iL^{2^)] (3) 

where luminosities are measured in erg s~^ and refer to the 
Ha line and the Spitzer 24-micron band respectively. 

We also estimate the to tal infrared luminosit y using all 
three MIPS bands following [Pale fc Heloul (|2002l ). 

i/TiR = wivL^{2A fiin)+W2i'L^{70 fim)+W3PL,, {160 jJ-m), (4) 

where wi = 1.559, 102 — 0.7686 and ws = 1.6381. The 
MIPS 160 /im band image is sampled without convolution 
since the resolution is larger than the HCN: 40" vs. 28", re- 
sulting in an additional uncertainty in Ltir of 15%. We also 
estimate the loc al pressure in the galaxy us ing the estima- 
tor presented bv lBlitz fc Rosolowskvl (|2006l ) and we report 
these values of P/k in Table [S] The total mass of the ob- 
served GMC as determined from the high-resolution map 
(R07) is given as Mco- The "Comparisons" section of the 
table highlights the most important results as discussed in 
O below. 



Figure 2. IRAM 30-m observations of HCN (1 0) and ^^CO 
(1 0) shown as black and gray lines respectively. The CO 
spectra have been scaled down by a factor of 100 for ease of 
comparison with the HCN data. The spectra have been Hanning 
smoothed and decimated to 2 km s^^ resolution. 



3 ANALYSIS AND RESULTS 

The minimal HCN emission from the galaxy was unex- 
pected given the predicted ratk)s of /co/.fHCN from other 
studies. For normal galaxies, ICao fc SolomonI (|2004h find 
Leo /Lhck ~ 25 w ith ~ 50% scatter. Fo r GMCs in the disk 
in the Milky Way, iHelfer fc Blit3 (|l997h find /co//hcn = 
40 ± 10 but the ratio decreases for clouds in the bulge of the 
galaxy to 12. For comparison, we find two sources (M33GMC 
91 and 26) with a dense gas ratio > 280 and the remaining 
two sources have a ratio below that of the Milky Way disk. 
We conclude that the HCN emission from M33 GMCs, rel- 
ative to their CO, is a factor of 2 to 7 lower than similar 
clouds in the disk of the Milky Way and other galaxies. 

The lack of HCN emission becomes even more remark- 
able when the results are compared with expectations based 
on the measu r ed sta r formation rate of M33. Previous work 
bv lWuetal.1 iooB^ suggests that the rate of ongoing star 
formation is linearly related to the luminosity of HCN (1 — >■ 
0): 



(Mq yr" 



1.4 X 10"''Lhcn (K km s" 



-1 2n 
pc ) 



(5) 



We estimate the amount of HCN emission that would 
be predicted based on the apparent star formation at each of 
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Table 2. M33 GMCs observed 



Property 


M33GMC 91 


M33GMC 26 


M33GMC 1 


M33GMC 76 


0:2000 


ni no 9 


ni "i*^ s 
ui 00 00.0 


ni f^9 A 


ni "iA 10 7 
Ul 0^ lU. / 


52000 


1 Qfi ziQ nf^ 


1 Qn A'\ 09 


T^ou oy 10 


iQn Qfi 1 


^LSR s ^) 


— z^y .u 


997 '\ 


1 fiQ n 
— luy .u 


1 f^Q 9 




246 


217 


244 


213 


a (mK) 


7.1 


3.1 


5.6 


5.5 


ho (K km s-i) 


21.6 


6.9 


7.2 


6.5 


/21 (K km s-i) 


25.2 


7.0 


15.5 


9.6 


/i3 (K km s-i) 


1.4 


0.3 


1.8 


1.6 


^HCN (K km s~^) 


< 0.053 


< 0.027 


~ 0.12 


~ 0.059 


LlQ (103 K km pc^)" 


150 


53 


47 


59 


Lhcn (10^ K km pc^)'' 


< 0.75 


< 0.40 


~ 1.90 


~ 0.94 


Ancillary Data 


-flO,R07 (K km s ^) 


29.0 


12 


11 


11 


.flO.convol (K km s~^) 


13.1 


5.0 


4.4 


5.2 




1.4 


2.9 


34 


10 


1/24^^^(24) (103 Lq)" 


61 


78 


470 


250 


Ltir (106 Ler 


1.4 


1.5 


6.0 


3.0 


SsFR (Mq Gyr-l pc-2) 


8 


12 


110 


43 


P/k (10"' K cm-3) 


1.6 


1.6 


2.4 


1.2 


Mco (105 Mq) 


11.1 


6.3 


3.1 


3.5 


Comparisons 


-^HCN (K km s~i) 


< 0.053 


< 0.027 


~ 0.12 


- 0.059 


-fnCN, pred (K km s"^) 


0.060 


0.090 


0.79 


0.310 


-flO,convol/-fHCN 


> 280 


> 190 


~ 37 


~ 91 


Ltir/ Lhcn'' 


> 2000 


> 4200 


~ 3300 


~ 3300 




7.2 


22 


92 


38 



° For conversion between intrinsic properties and quantities derived from surface brightness, note 
that the 28" IRAM beam subtends 1.5 X 10* pc^ at the 840 kpc distance of M33 llFreedman et al.l 

[iooj). 

Lq/( K km s-i pc2) 



our pointing locations in M33. Using the star formation rate 
estimates presented in Table (2] we use Eq. [5]to predict the 
intensity of HCN emission expected from the ongoing star 
formation (/ncN.pred in Table[2]). We find the predicted HCN 
emission is a factor of 1.1 to 6.6 higher than the observed 
values or upper limits. The discrepancy may result, in part, 
from the inclusion of Ha in the star formation indicator vs. 
the infrared-only data used in Equation (5) 

The "CO (2 1) and ^=^00(1 -> 0) data ofTer 
some insight into the conditions of the molecular cloud 
targets. These line ratio data can also compared to the 
Ao,R07 data. The difficulty in evaluating these results is 
the difference in beam sizes between the ^^CO(2 1) line 
(11") and the "C0(1 -5> 0) line (14"). If the emission 
fills all beams, then the line ratio can be m odelled using 
LVG analysis as was done for M33 GMCs in I Wilson et al.l 
( 1997h. We attempt an LV G analysis using the RADEX code 
( van der Tak et al]|2007t ) . We take a typical column density 
oiNco ~ 10^'^ '' cm-2 (assuming iV(CO)/Af(H2 = 2x 10"^) 
and line FWHM of ~ 8 km s~^ representative of the sample. 
We a dopt iV("CO)/iV("CO)=70 following Iwilson et all 
l|l997h though they argue that LVG results are relatively 
insensitive to this ratio. We calculate an LVG model for 
two extreme conditions: (1) the emission is uniform within 
all the beams so the source-beam coupling is unity or (2) 
the emission is a point source within the beam so the in- 



tensities must be scaled by the ratio of the beam areas 
for comparing the lines. We model the expected line ratios 
"CO(2 ^ l)/"CO(l ^ 0) and ^='CO(2 ^ 1)/"C0(2 ^ 1) 
for a grid of density (n = 10^'^ — >■ lO^'^ cm"'') and temper- 
ature conditions {Tk = 5 K — >■ 100 K). We find solutions 
where the physical conditions produce a consistent set of 
line ratios and label those points with the M33GMC num- 
ber in Figure [3] The models have no free parameters and 
thus are able to replicate line ratios exactly, provided the 
results are sampled on the grid. Observations of more lines 
would enable a statistically robust determination of cloud 
properties. While the permitted range of conditions is broad 
depending on assumptions regardings source structure, we 
find no cloud that is inconsistent with all clouds being typ- 
ical GMCs (i.e., n ~ 10=* cm'^ and T ~ 15 K). While there 
is some concern that the conditions in the clouds may not be 
well suited for HCN (1 — >■ 0) emission, the line ratio studies 
are consistent with these objects being normal GMCs. 



4 DISCUSSION 

The minimal HCN emission from M33 is striking, especially 
compared with the active star formation from the galaxy 
as a whole. In particular, M33GMCs 91 and 26 highlight 
the deficit of HCN emission relative to CO when compared 
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Figure 3. Large Velocity Gradient analyses of the pointings in 
the observation using the RADEX code applied to the ^'^CO{2 -i> 
1) and (1 0) lines and the ^^CO{2 1) line. The top and 
bottom panels indicate the result assuming that all the emission 
from the GMC targets is either a point source (top) or beam- 
iilling (bottom). The lines are shown to indicate conditions pro- 
ducing a line ratio. Intersections between pairs of lines show where 
conditions produce both line ratios, which are labelled with the 
corresponding M33GMC number. Points occur at the edge of the 
model grid if the solution is outside the sample box. 



to Milky Way clouds throughout the disk. We note that 
M33GMC 91 is the most massive cloud in M33. M33GMCs 
1 and 76 show marginal detections consistent with the outer 
bounds of what is expected from Milky Way and other local 
galaxies. In this section, we explore possible explanations for 
the high CO/HCN ratios observed. 

The best comparable study of CO-to-HCN ratios in 



Figure 4. CO-to-HCN ratios for low-z galaxies ijCao fc Solomon! 
I2OO4I . , grey) and for M33 pointings (black). The distribution 
shows that the HCN emission is anomalously weak compared to 
CO when viewed in the context of galaxies. The highest value 
of the CO-to-HCN emission is associated with the most massive 
GMC in M33. 

Local Group galaxies is the work of lBrouillet et al ] (|2005l . 
BOS), who used the IRAM 30-m telescope to observe 16 
lines of sight in M31. Their work found a radial decrease 
in 7co /Ihcn ranging from 30 to 125 at the outer limit (see 
Figure [Jjl . BOS interpreted the radial variation as a reflec- 
tion of a decrease in dense gas content of molecular clouds 
at large r galactocentri c radiu s, consistent with the decrease 
seen in iHelfer fc Blit j l| 19971 ). However, both M31 and the 
Milky Way show relatively normal star formation efficien- 
cies for their molecular gas content, in contrast with M33. 
Further, all of the B05 m easurements are i n regions with 
^coZ-'^HCN and metallicity (jBlair et al.lll983 ) similar to the 
Milky Way measurements. 

4.1 Excitation Conditions 

In this work, we have focused our analysis on the (1 — >■ 0) 
lines of ^^CO and HCN. Following a significant number of 
authors, we have taken the CO emission to trace the "low- 
density" emission in a molecular cloud (n ~ 10^ cm""^), 
whereas the HCN line traces high density gas within the 
cloud. The latter statement is largely motivated by the high 
critical density of the transition. However, the critical den- 
sity is only a rough indicator of the excitation conditions of 
the gas. High temperatures (T > 50 K) and significant ra- 
diative pumping can bring sharp departures from this "com- 
mon sense." Without a full, multi-line analysis of the species, 
a firm interpretation is impossible {^^. However, the anal- 
ysis of the physical conditions in M33 clouds suggest that 
these sys tems are a nalogous to GMCs found in the Milky 
Way (.Wilson et al][l997 ; Heifer fc Blit3ll997l ). rather than 
active galactic nuclei where substantial concerns about the 
interpretation of HCN have been raised. Hence, relying on 
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the (1 — !> 0) lines of CO and HCN as tracers of low- and 
high-density gas respectively appears justified. 

4.2 Metallicity and Chemical Abundances 

The M33 system is sub stantially more metal poor 
l|Rosolowskv fc SimonI [20081 ) with all but the central mea- 
surement in the present study occurring at lower metallicity 
than any point in the M31 stud y. The metallicity m ay ap- 
proach solar in the inner 1 kpc ( Magrini et al.ll2007l ) which 
harbours M33GMC1. The similarities in pressure (see be- 
low) and the presence of significant amounts of star forma- 
tion in the clouds suggests that there should be sufficient 
dense gas to excite HCN emission if it is present. The lower 
overall metallicity is also related to depletion of secondary 
(N) elements relative to primary elements (C,0), so HCO-I- 
may be a more viable tra c er of dense gas in low metallic- 
ity svstems. [Magrini et al] |2007) find a [N/O] ratio 0.3 dex 
lower than solar abundances across the entire galaxy. This 
depletion may translate directly to lower HCN/H2 abun- 
dances. 

To evaluate the effects of lower metallicity on the rel- 
ative abundance of chemical species, we i ntegrated a set 
of ch emical models using the Nahoon code jWakelam et al.l 
|2004 ). We considered a set of eight models spanning a range 
from log[(O/H)/(O/H)0] = 0.0 to -0.4 in uniform logarith- 
mic steps and log[(N/O)/(N/O)0] = 0.0 to -0.3 in uniform 
logarithmic steps. Solar abundances were set at the Nahoon 
defaults: He/H = 0.14, C/H = 7.3x 10"^ , N/H = 2.14x 10-^ 
0/H = 1.76 X 10"'', Cl/H = 4 X 10"^ F/H = 6.68 x 10"^ 
Fe/H = 3 X lO-'-*, Mg/H = 7 x 10"^ Na/H = 2 x 10-^ P/H 
= 3 X 10-^ S/H = 8 X 10"*, Si/H = 8 x 10"^ We scaled 
the solar abundances of a-process elements to track the O 
abundance (C, Mg, S, Si) and all other elements to track N 
(CI, F, P, Na). Grain abundance was depleted by the same 
factor as O. We use standard initial conditions for chemi- 
cal models, including nH2 = 10^ cm"'^ and T = 10 K. The 
models use the 2005 OSU Astrochemistry Rate compilation 
distributed with the Nahoon package. The rates track 4423 
reactions among 452 species. 

We compared the abundance of several chemical species 
at t = 10^ years (the typical lifetime of a CMC) to evaluate 
the effects of depletion and low metallicity on the forma- 
tion of various species, specifically HCN and CO. We find 
the CO abundance is a reasonable value for the solar abun- 
dance model (CO/H2 = 1.4 x 10~*) and decreases roughly 
linearly with the decreasing metallicity. Contrary to expec- 
tations, the HCN abundance is nearly constant with decreas- 
ing metallicity (HCN/H2 = (1.4±0.1) x 10"^), as a resuh the 
CO/HCN ratio actually decreases by a factor of 2.5 going 
from solar to M33 abundances. These models are necessar- 
ily coarse approximations to the actual chemistry present in 
clouds. However, they do not indicate significant increase of 
the CO/HCN abundance with a decreasing N/O ratio. The 
ratio HCO-I-/HCN is also roughly constant with metallicity. 
While not conclusive, these models are a reminder that the 
abundances of individual chemical species need not track the 
abundances of elements of which they are composed. Hence, 
metallicity chan ges do not app ear to explain the dearth of 
HCN in M33. iHilv-Blant et all |2010.) have shown that the 
HCN ratio is insensitive to C/O variations in the gas, similar 
to the model's N/O and O/H robustness. A full treatment, 



including the effects of radiation ( as was done for CO in 
iBell et al.1 [2OO6I : IClover et al] 120101 ) . is needed to validate 
this initial exploration. 

Some authors have raised concerns about the utility 
of HCN as a star formation tracer, noting that AGN ac- 
tivity could decrease the CO/HCN ratio such that there 
would appear to be more dense gas than expected (e.g. 
iMeiierink et al.ll2007l ). We note that this is not a concern 
here because such effects invoke x-ray irradiation by an AGN 
and such effects should make HCN easier rather than more 
difficult to detect. 

Studies of HCN formation in photodissociatio n regions 
([Young Owl et al.1 12OO0I : iBoger fc Sternberg! |2005| ) suggest 
that HCN is readily photodissociated in translucent regions. 
Their work shows that HCN does not reach significant con- 
centrations unless shielded by 4 to 8 magnitudes of visual 
extinction, depending on the local conditions. Given that 
CO forms more readily [Av ^ 2), translucent clouds will 
show a high CO/HCN ratio. If the subsolar metallicity in 
M33 predicts for lower dust-to-gas ratios, then the fraction 
of the volume of the molecular cloud at high Ay may be sig- 
nificantly reduced, increasing the CO/HCN ratio. Identical 
reasoning is forwarded to explain the low brightness of CO 
in the SMC and other low metallicity systems (e.g., Ilsraell 
ll997l : lLerov et al.ll2009l . among others) . Given the ample ev- 
idence for such an effect in low metallicity systems, it may 
hold that HCN will cease to be a good tracer of star for- 
mation at modest metallicities (e.g., log([0/H]) -I- 12 ~ 8.4) 
because of dissociation effects. Further studies of the dust 
content of M33 are ongoing and should clarify the impor- 
tance of photodi ssociation at regul ating the CO/HCN ratio 
on CMC scales ijBraine et al.ll2010l ). 

4.3 ISM Structure 

Changes in the structure of the ISM are important at setting 
the CO/HCN ratio in starburst systems, where the dense 
gas o ccupies a large fraction of the ISM (|Gao fc Solomon! 
!2004 ). This change reduces the CO/HCN ratio. It is con- 
ceivable that a similar effect is at work here. However, the 
molecular cloud population in M 3 3 has been studied ex- 
tensi vely (|Wilson fc ScovilQ !l990l : !Rosolowskv et aL I !200l 
!2007h . reaching the conclusion that the individual clouds 
have similar macroscopic properties as those found in the 
Milky Way and M31, including turbulent line widths and 
average densities. 

The properties of these clouds are linked to t he galac- 
tic environment via pressure. !Helfer fc Blitl (|l997! ) demon- 
strate that HCN-to-CO ratio correlates with the disk pres- 
sure in the ISM at the galactocentric radius of the ob- 
jects observed: Ihcn/J'co oc pO i9±o o-i Higher pressures in 
the galactic environment would cause GMCs to have cor- 
respondingly higher pressure to maintain their distinction 
from the remainder of the ISM. As seen in Table ^ the 
disk pressure values are comparable to the Milky Way disk 
va lues, implying these GMCs do not fall on the scalings 
of !Helfer fc Blitg (|l997! ). Following their scaling, we predict 
^co/-^HCN ~ 30 — 50 in these clouds which is only consis- 
tent with a marginal detection from M33GMC 1. Thus, it 
seems unlikely that the dense gas content of these clouds 
would be significantly lower than those in the other local 
group disk systems. This is because the dense gas fraction 
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of clouds has a tight link betwe en the average density and 
amount of turbulent driving (e .g. JPadoan fc Nordlundl2002l : 
iKrumholz fc Thompsonll2007l ). 



4.4 Infrared Luminosity and Evolution 

The motivation for studying HCN in M33 comes from the 
intimate link between star formation (as implied by infrared 
emission) and HCN emission in galaxies. In Figure [3 we 
consider the infrared and HCN emission in the context of 
other studies. 

The targets in M33 are comparable to high luminos- 
ity star forming regions in the Milky Way. While they ap- 
pear significantly above the mean ratio defined by galax- 
ies (I/ir/I/hcn = 800), the lines-of-sight appear consistent 
with the upper envelope of values derived for the Milky 
Way regions. We note that t he observed clouds have sub- 
stantial masses (> lO'^ Mm. iRosolowskv et al.1 [ioOTi ) and 
their infrared luminosities are well above the threshold above 
which the HCN-to-IR rat io saturates at its constant ratio 
(Figure [51 IWu et al.lboicl '). Specifically, M33GMC 1 and 76 
have bright, discrete Hll regions located within the IRAM 
beam, in additi on to being strong CO emitters. Although 
IWu et al.l l|2010l ) have argued for a constant Lir/Lhcn ra- 
tio for LiR > 10*'^ L0, the addition of the M33 points to the 
upper envelope of the trend suggest a continuation to higher 
ratios at larger luminosities. Hence, the remarkably consis- 
tent values for the IR-to-HCN ratio seen in galaxies may 
actually derive from consistent sampling from a population 
with a changing ratio. 

The relationship between star formation and dense 
molec ular gas emission is not universal. Indeed, ICao et al] 
(1200711 found that four high redshift galaxies {z = 1.0 2.8) 
show LncN/Af* that is a factor of 2 lower than the result es- 
tablished in the local population and several non-detections 
of HCN suggest a cosmological evolution of the star forma- 
tion efficiency of dense gas. While the GMCs in M33 do not 
have star formation efficiencies as high as observed in some 
high redshift systems, the results are consistent with the star 
formation efficiency of HCN-emitting gas being comparable 
to or larger than that observed in high redshift systems. 

This study observes i ndividual GMCs whereas the other 
studi es focus on galaxies l|Gao fc Solomonl l2004l: ICao et akl 
[2003) or individual star forming regions jWu et al.ll2005h . 
Thus, the evolution of GMCs via cloud evolution and star 
formation processes may cause scatter in the interpretation 
of star formation efficiency. We have selected GMCs such 
that I/ir/I/co spans a large range. While not conclusive, it 
is suggestive that the two marginal detections showing the 
lowest IcoZ-'^HCN ratios are also the most infrared luminous. 
The lack of dense gas in M33GMCs 91 and 26 may indicate 
these are younger, less evolved systems. The factor-of-ten 
range in star formation rates derived from the IR -I- Ha 
emission suggests that variations in evolutionary stages may 
indeed be responsible for the low levels of HCN emission seen 
from M33GMCs 91 and 26. Thus, a possible explanation for 
the highest values of CO-to-HCN may be that these clouds 
are young and have not formed significant fractions of high 
density gas. 



5 CONCLUSIONS 

In summary, we have searched for and failed to find sig- 
nificant HCN emission from GMCs in M33. The HCN ob- 
servations are sufficiently sensitive that the HCN should be 
easily detected (> 20crrms) if the emission were present at 
the levels expected based off either CO or star formation 
observations in our Galaxy and other s (Hclfcr & Blitz 1993i 
iGao fc Solomonll2004l : IWu et aLllioOsI ) . Two of the four lines 
of sight that we observed show a marginal detection of HCN 
emission and these two targets are bright infrared sources 
in the galaxy. Using naive chemical models, we explored the 
effects of the low metallicity and N/O depletion in M33 on 
the CO/HCN ratio. Despite depletion of N/O, the CO/HCN 
ratio is predicted to be lower in M33 than in the Milky Way. 
Since HCN only forms at high extinction relative to C/O, 
the low metallicity of the galaxy may result in reduced dust 
shielding and thus more molecular gas exposed to dissoci- 
ating radiation. Even though the ratio of -Lir/Lhcn is a 
factor of six larger than the mean for galaxies, it is nonethe- 
less consistent with comparable systems in the Milky Way. 
The M33 targets appear most distinct in their CO-to-HCN 
ratios rather than in their infrared luminosities. This sug- 
gests that the low levels of HCN emission likely reflect a 
smaller fraction of dense, shielded gas in these clouds rather 
than a systematic effect. Having considered chemical and 
structural variations, the observed discrepancies seem most 
likely to be the result of either reduced extinction or cloud 
evolution effects. A larger sample of GMC-based observa- 
tions and a clearer understanding of the dust properties in 
M33 are needed to clarify these results. 
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